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Abstract. A quantitative understanding of the weak nuclear response is a prerequisite for 
the computer simulations of astrophysical phenomena like supernovae explosions and neutron 
star cooling. In order to reduce the systematic uncertainties associated with the simulations, 
a consistent framework, able to take into account dynamical correlation effects, is needed to 
compute neutrino-nucleon and neutrino-nucleus reaction rates. In this paper we describe the 
many-body theory of the weak nuclear response at low energy regime. We show how to include 
both short and long correlations effects in a consistent fashion. 



1. Introduction 

The description of neutrino interactions at low momentum transfer (E u ~ 10 MeV) with 
nuclei, and nuclear matter in general, is relevant to the study of many different problems, 
from supernovas explosions [1] to neutron star cooling [2]. 

The systematic uncertainty associated with computer simulations depends heavily on the 
values of the neutrino-nucleon and neutrino- nucleus reaction rates used as inputs. 

Nuclear many body theory provides a scheme allowing for a consistent treatment of neutrino- 
nucleus interactions. Within this approach, nuclear dynamics is described by a phenomenological 
hamiltonian, whose structure is completely determined by the available data on two- and three- 
nucleon systems, and both short and long range dynamical correlations are taken into account. 

Over the past decade, the formalism based on correlated wave functions, originally proposed 
to describe quantum liquids [3], has been employed to carry out highly accurate calculations of 
the binding energies of both nuclei and nuclear matter, using either the Monte Carlo method 
[U \E\ E] or the cluster expansion formalism and the Fermi Hypernetted Chain integral equations 
[ZlEli- 

A different approach, recently proposed in Refs. |10|lll] exploits the correlated wave functions 
to construct an effective interaction suitable for use in standard perturbation theory. This scheme 
has been employed to obtain a variety of nuclear matter properties, including the neutrino mean 
free path [10J and the transport coefficients [111 I12j . 

In this work we describe the application of the formalism based on correlated wave functions 
and the effective interaction to the calculation of the weak response of uniform nuclear matter. 



2. Many-body theory of the weak nuclear response 

2.1. Neutrino-nucleus cross section 

The differential cross section of the process 

v t + A -» I + X , (1) 

in which a neutrino carrying initial four- momentum k = (E u ,\s) interacts with a nuclear target, 
producing a lepton in a state of four-momentum k' = (E^,k ; ), the target final state being 
undetected, can be written in Born approximation as (see, e.g., Ref. [16] ) 
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where G = Gf cos 9c, Gf and #<7 being Fermi's coupling constant and Cabibbo's angle. The 
leptonic tensor, that can be written, neglecting all lepton mass, as 



kjxk v + k u k /1 9[n>{kk) i^^ua/3^ , (3) 



is completely determined by lepton kinematics, whereas the nuclear tensor W^ v contains all the 
information on target structure. Its definition involves the initial and final hadronic states |0) 
and \X), carrying four- momenta po and px, respectively, as well as the nuclear weak current 
operator J^: 

W>» = J2(0\J tl \X)(X\r\0)6^{p + q-p x ) , (4) 
x 

where the sum includes all hadronic final states. 

In the low momentum transfer regime, we expect the non relativistic approximation to be 
applicable. Within this approach, the initial and final states can be obtained from nuclear many- 
body theory (NMBT), while the weak current entering the definition of the hadronic tensor [1] 
is expanded in powers of |q|/m. At leading order, the resulting response can be written in the 
simple form 

S M = ^E( l°ql n )( n l°ql )^ + E ° - E n) ■ (5) 

n 

where, in the case of charged current interactions, O q is the operator corresponding to Fermi or 
Gamow- Teller transitions, whose expressions in the coordinate space are: 



0[( q )=g v S(r i -r' i ) e^r+ 
OF T (q)=^(r i -r^)e^<x l T+, (6) 



where r^ specifies the position of the i— th particle. 



2.2. Description of the initial and final states 

Understanding the properties of matter at densities comparable to the central density of atomic 
nuclei is made difficult by both the complexity of the interactions and the approximations implied 
in any theoretical description of quantum mechanical many-particle systems. 

The main problem associated with the use of the nuclear potential models in a many-body 
calculation lies in the strong repulsive core of the NN force, which cannot be handled within 
standard perturbation theory. 



Within NMBT, a nuclear system is seen as a collection of point-like protons and neutrons 
whose dynamics are described by the hamiltonian 



# = £i» + £^j) + ... , (7) 

i j>i 

where t(i) and v(ij) denote the kinetic energy operator and the bare NN potential, respectively, 
while the ellipses refer to the presence of additional many-body interactions. 

Carrying out perturbation theory in the basis provided by the eigenstates of the 
noninteracting system requires a renormalization of the NN potential. This is the foundation of 
the widely employed approach developed by Bruckner, Bethe and Goldstone, in which v(ij) is 
replaced by the well-behaved G-matrix, describing NN scattering in the nuclear medium (see, e.g. 
Ref.|13j). Alternatively, the many-body Schrodinger equation, with the hamiltonian of Eq. ([7]), 
can be solved using either the variational method or stochastic techniques. These approaches 
have been successfully applied to the study of both light nuclei [3J and uniform neutron and 
nuclear matter [71 151 IS} [T7] . 

Our work has been carried out using the scheme of correlated basis function (CBF) theory 
in which nonperturbative effects due to the short-range repulsion are embodied in the basis 
functions. 

The correlated states of nuclear matter are obtained from the Fermi gas (FG) states \ufg) 
through the transformation [3l [18] 

.F|bfg) 

W = (n FG |FtF|n FG )V2 ' ^ 

In the above equation, |tt-fg) is a determinant of single particle states describing iV noninteracting 
nucleons. The operator F, embodying the correlation structure induced by the NN interaction, 
is written in the form 

N 

F(1,...,N)=S J] f i:i , (9) 
j>i=i 

where S is the symmetrization operator which takes care of the fact that, in general, 

[fijJi^^O. (10) 

The structure of the two-body correlation functions must reflect the complexity of the 
NN potential. The shape of the radial functions f n (rij) is determined through functional 
minimization of the expectation value of the nuclear hamiltonian in the correlated ground state 

EX = (0\H\0) . (11) 

As an example, FigQ] shows the radial dependence of the potentials ) and correlation functions 
acting in the spin-isospin channels S = and T = and S = and T = 1. 

It has to be pointed out that the correlation operator of Eq.([9|) is defined such that, if any 
subset of the particles, say i\, . . . i p , is removed far from the remaining i p +\, . . . *jv> it factorizes 
according to 

F(1,...,N) ^ F p (i u ...i p )F N _ p (i p+1 ,...i N ) . (12) 

The above property is the basis of the cluster expansion formalism, that allows one to write 
the matrix element of a many-body operator between correlated states as a sum, whose terms 
correspond to contributions arising from isolated subsystems (clusters) involving an increasing 
number of particles. 




Figure 1. Interaction potentials (upper panel) and correlation functions (lower panel) acting 
in the spin-isospin channels S = and T = (solid lines) and S = and T = 1 (dashed lines). 
The potential is the Argonne v' 8 [15] and the correlation functions correspond to nuclear matter 
at equilibrium density. 



3. Correlation effects in the transition matrix elements 

Using correlated states implies severe difficulties in the explicit calculation of the weak matrix 
element. In the FG model, the nuclear response is non vanishing only when the final nuclear 
state differs from the initial state for the presence of a particle excited outside the Fermi sea 
and a hole in the Fermi sea. In the presence of correlations, which can induce virtual nucleon- 
nucleon scattering processes leading to excitation of nucleons to states outside the Fermi sea, 
more complex scenarios must also be considered. For example, if the initial state has a two 
particle-two hole component, the final state can be a three particle-three hole state or, if the 
probe interacts with an excited nucleon, a two particle-two hole state. 

In the following we will consider only the dominant transition, between the correlated ground 
state and a correlated one particle-one hole (ph) state. The corresponding weak matrix element 
can be written 

M vk = W F ? F '°> , , (13) 

{ph\FiF\ph)*(0\FiF\Q}2 

where F is the correlation operator defined in Eq.(|9|). Here the kets |0) and \ph) correspond to 
the ground and one particle-one hole Fermi Gas states, respectively, and O = J2iOi, Oi being 
the Fermi or the Gamow- Teller transition operator (see Eqs.([6])). 

Note the the = fij — 1 is short ranged, and therefore its matrix elements are small. At 
two-body level, the cluster expansion of Eq. (fT3j) yields [20] 

(ph\F*OF\0) ~ (ph\(l + J2 9^)0(1 +Y, 9 ij)\°) ■ ( 14 ) 

j>i j>i 

The above equation suggests the definition of an effective operator OyJ, acting on Fermi Gas 
states. From 

(ph\F^OF\0) = {ph\O eff \0). (15) 



it follows that, at the two-body cluster level (compare to Eq. (ll4|) ) 
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where {A, B} = AB + BA and N denotes the number of particles. Note that the O e ^ is a 
two-body operator, as it includes screening effects arising from nucleon-nucleon correlations. 

We have carried out the calculation of M p h on a cubic lattice, for a discrete set of states 
hj, satisfying the conditions hj < kf and h{ + q > kf, kf being the Fermi momentum. Figure 
[2] shows a comparison between the correlated and the Fermi Gas matrix element for a Fermi 
transition. It has to be noticed that the effects of correlations is enhanced in the calculation of 
the response, see Eq. ([!)]), whose definition involves the square of the correlated matrix element. 
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Figure 2. Fermi transition matrix element at |q| = 0.3 fm 1 as a function of the magnitude of hole 
momentum |h|. The dashed horizontal line corresponds to the result of the FG model. 



4. The effective interaction 

At lowest order of CBF, the effective interaction V e R is defined by 

(H) = (0 FG \F^HF\0 FG ) = (0 FG \T + V eS \0 FG ) . (17) 

As the above equation suggests, the approach based on the effective interaction allows one to 
obtain any nuclear matter observables using perturbation theory in the FG basis. However, as 
discussed in the previous Section, the calculation of the hamiltonian expectation value in the 
correlated ground state, needed to extract V e g from Eq. (fTTl) . involves severe difficulties. In this 
work we follow the procedure developed in Refs. [HIHO], whose authors derived the expectation 
value of the effective interaction by carrying out a cluster expansion of the rhs of Eq. (|17p . and 
keeping only the two-body cluster contribution. The resulting expression can be written 

{0 FG \V cS \0 FG ) = Y^(ij\Ve&(12)\ij)a = X>j|/l2 

i<j i<j 

where the laplacian operates on the relative coordinate and the suffix a denotes that the ket 
\ij) a is anty-simmetryzed. 

In the left panel of Fig|3] the central components of the effective interaction obtained from 
the v' 8 potential of Ref.[T5] at equilibrium density is compared to the corresponding component 



— (V 2 /l2)+^(12)/l2 
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Figure 3. Left panel: Comparison between the central component of the bare Argonne fg[T5] potential 
(dashed lines) and the effective potential defined by Eq. (fT7|) (solid lines), calculated at nuclear matter 
equilibrium density. Right panel: momentum dependence of the single particle energies obtained from 
the CBF effective interaction in the Hartree-Fock approximation (solid line). The diamonds show the 
results of Ref.|19j. For reference, the energies of the Fermi Gas model are also shown (red line): 

of the bare interaction. It clearly appears that screening effects due to NN correlations lead to 
a significant quenching. 

Using the effective interaction described above, the single particle energies ek in nuclear 
matter can be easily computed in Hartree-Fock approximation. The risulting expression is: 



The right panel of Fig. [3] shows the energy spectrum for symmetric nuclear matter at 
equilibrium density. 

5. Correlation effects on the response. 

5.1. Short range correlations. 

We have calculated the nuclear response for a Fermi transition (Eqs.© and ©)• The inclusion 
of interaction leads to sizable modifications of the FG response. Correlation effects in the 
transition matrix elements, taken into account through the use of the effective operator, produce 
a quenching of ~ 15 %. This feature is apparent in Fig. [H where the FG response, represented by 
the green line, is compared to that obtained using the effective operator in the calculation of M v h 
and the FG single particle spectrum, represented by the red line. An even larger modification is 
produced by interaction effects on the single particle energies. Replacing the FG single particle 
energies with the HF energies (Eq ll9p leads to a sizable broadening of u region corresponding 
to non- vanishing response as shown by the black line in Fig. [H 

5. 2. Long range correlations: the Tamm-Dancoff Approximation ( TDA ) 

In the previous section we have discussed the nuclear response in the correlated Hartree-Fock 
(HF) approximation, in which the bare Fermi transition operator is replaced by the effective 
operator of Eq. (|15p and the final state is assumed to be a one particle-one hole state. 

It is important to realize that the FG one particle-one hole states, while being eigenstates of 
the HF hamiltonian, defined as 
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Figure 4. Nuclear response for Fermi transition, for q = .3 fm _1 (left), q = 1.8 fm _1 (center), 
and q = 3 fm^ 1 (right). Black line is full calculation. Confrontation with Fermi Gas (green 
line) and effective operator with Fermi Gas energy spectrum (red line) are shown. 

with ej given by Eq, (|19|) . are not eigenstates of the full nuclear hamiltonian. As a consequence, 
there is a residual interaction V res , which can be identified with the effective interaction defined 
in Section HI that can induce transitions between different one particle-one hole states, as long 
as their total momentum q, spin and isospin are conserved. 

In order to include the effects of these transitions, we use the TDA, which amounts to 
expanding the final state in the basis of one particle-one hole states according to|21j 



|/) = |q, TSM) = ]T cJ SM \h h Pl = hi + q, TSM) , (21) 

i 

where S and T denote the total spin and isospin of the particle hole pair and M is the spin 
projection. 

At fixed q, the excitation energy of the state (|21|) . oj? , as well as the coefficients c[ SM , are 
obtained solving the eigenvalue equation 

H\f) = (H HF + V res )\f)= (» • X>e//fa')l I/) = (^0 + ^)1/) , (22) 
where Eq is the ground state energy. In the TDA, the response can be written as 



T,M,S n 



J2(c T n SM h(hi, p h TSM\O eff (q)\0) 
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TSM 



5(u-^), (23) 



where {c^ SM )i denotes the z-th component of the eigenvector belonging to the eigenvalue ■ 
We have performed the diagonalization of the Hamiltonian operator over a basis of ~ 3000 
one particle-one hole states for any spin-isospin channel. The appearance of an eigenvalue lying 
well outside the particle hole continuum, corresponding to a collective excitation, reminiscent 
to the plasmon mode of the electron gas, is clearly visible in Figj5j where the calculation of the 
response for a Fermi transition within TDA is compared with the response of the correlated HF 
approximation. The TDA response exhibits a sharp isolated peak corresponding to the collective 
mode, lying ~ 4 MeV above the upper limit of the particle hole continuum. 



6. Conclusions 

We have carried out calculations of the charged current weak response of nuclear matter in the 
low momentum transfer regime. The quantitative understanding of this quantity is required in 
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Figure 5. Nuclear matter response in the TD approximation at |q| = 0.3 fm 1 (solid line). The dashed 
lines show the results of the correlated HF approximation. 

many problems of physics, ranging from simulations of supernovas explosions to neutron star 
cooling. 

The calculation has been performed using a many-body approach based on a realistic nuclear 
hamiltonian, including two- and three-nucleon interactions, yielding a good description of 
the properties of both the two-nucleon systems and uniform nuclear matter. The response 
associated with Fermi transitions at low momentum transfer has been calculated from an effective 
interaction, derived using the formalism of correlated basis functions and the cluster expansion 
technique. Our work improves upon existing effective interaction models in that it includes the 
effects of many-nucleon forces, which become sizable, at high density. The responses calculated 
within the correlated HF approximation show that the inclusion of short range correlations leads 
to a significant quenching of the transition matrix elements and shifts the strength towards larger 
values of the energy transfer, u>, for all values of |q|. At low momentum transfer, long range 
correlations have been taken into account within the TDA. The excitation of the coherent state 
can be clearly seen in our results at |q| = 0.3 fm -1 . 
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